Short-time inertial response of viscoelastic fluids measured with Brownian motion and 

with active probes 
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We have directly observed short-time stress propagation in viscoelastic fluids using two optically 
trapped particles and a fast interferometric particle-tracking technique. We have done this both 
by recording correlations in the thermal motion of the particles and by measuring the response 
of one particle to the actively oscillated second particle. Both methods detect the vortex-like flow 
patterns associated with stress propagation in fluids. This inertial vortex flow propagates diffusively 
for simple liquids, while for viscoelastic solutions the pattern spreads super-diffusively, dependent 
on the shear modulus of the medium. 

PACS numbers: 83.60.Bc,66.20.+d,82.70.-y,83.50.-v 



I. INTRODUCTION 



Motion in simple liquids at small scales is usually char- 
acterized by low Reynolds numbers, in which the re- 
sponse of a liquid to a force applied at one point is Stokes- 
like — decaying with distance r as 1/r away from the ori- 
gin of the disturbance [3, H, [H . Here, fluid inertia can 
be neglected, and the force is effectively felt instanta- 
neously everywhere within the medium. In practice, this 
is a good approximation, for instance, in water at the 
colloidal scale up to micrometers on times scales larger 
than a microsecond. At short times or high frequencies, 
however, fluid inertia limits the range of stress propaga- 
tion. Any instantaneous disturbance must be confined 
to a small region for short times. If the medium is also 
incompressible, then this naturally gives rise to vorticity 
and backflow. In simple liquids, stress then propagates 
through the diffusive spreading of this vortex. Although 
this basic physical picture has been known theoretically 
for simple liquids since the work of Oseen in 1927 [4[ 
and has been shown in computer simulations since the 
1960s [f|, experimental observation of these effects have 
been largely indirect, for instance, in the form of short- 
time corrections to Brownian motion [f|. Direct exper- 
imental observation has only recently been possible be- 
cause of the high temporal and spatial resolution required 

0- 

The finite time it takes for vorticity to propagate leads 
to persistence of fluid motion that manifests itself in al- 
gebraic decay of the auto-correlation function of the ve- 
locity of either a fluid element or a particle embedded in 
the fluid. This decay is slower than the naive expectation 
of exponentially-decaying correlations for a massive par- 
ticle experiencing viscous drag. Thus, the effect is known 



as the long time tail effect [6j, |8|, |9( , which characterizes 
the transition from ballistic to Brownian motion of par- 
ticles in simple liquids. This effect has been shown to 
be present even at the atomic level, e.g., from neutron 
scattering experiments on liquid sodium [10| . 

We have shown that the inter-particle correlations and 
response functions of two particles can be used to di- 
rectly resolve the flow pattern and dynamics of vortex 
propagation 0, El . This was done by measuring the cor- 
related thermal motion of two optically trapped spheri- 
cal particles using an interferometric technique [la [l3j 
with high temporal (up to 100 kHz) and spatial (sub- 
nanometer) resolution in both viscous and viscoelastic 
fluids. We were able to observe, for instance, the anti- 
correlation in the inter-particle fluctuations of thermal 
motion that is characteristic of the vortex propagation. 
The method is related to passive two-particle microrheol- 
ogy [3, [H, [H, E3, Gil EH which can be used to measure 
shear elastic moduli of viscoelastic materials. 

The inter-particle response functions a(u), with real 
(a') and imaginary (a") parts defined by a(u>) = a'(cu) + 
ia"(to) are obtained for motion parallel and perpen- 
dicular (a±) to the center line connecting the two parti- 
cles. In the passive approach, we directly measure the 
imaginary part of the response function from the ther- 
mal position fluctuations of the two particles via the 
fluctuation-dissipation theorem (FDT). The real part of 
the response function is then obtained from a Kramers- 
Kronig integral [2(| ■ 

Here, in order to directly measure both real (in-phase) 
and imaginary (out-of-phase) parts of the res pon se func- 
tion, we have developed an active method [2 lL [22l |. in 
which one optical trap drives oscillatory motion of one 
particle, while the response of a second particle is mea- 
sured at separation r. We present and compare de- 
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tailed experimental results of both passive and active 
approaches. We also present a theoretical derivation of 
the predicted response functions and corresponding alge- 
braic decay of the velocity autocorrelation functions for 
viscoelastic fluids. 

The outline of the paper is as follows. In section II, 
we present the theoretical analysis. In section III, we 
present the materials and methods of sample prepara- 
tion, as well as the experimental techniques for the pas- 
sive and active measurements of the response functions. 
In section IV we describe our methods of data analysis 
used for the results presented in section V. In the results 
section V, we first compare the data for simple liquids 
with the dynamic Oseen tensor, which demonstrates the 
diffusive propagation of the vortex flow. We then present 
our results and comparison with theory for viscoelastic 
solutions, including the evidence for superdiffusive stress 
propagation. Finally, we conclude with a discussion (sec- 
tion VI), also mentioning implications of our results for 
microrheology in general. 



dependence of the velocity field. For an oscillatory distur- 
bance at frequency u>, the propagation of vorticity defines 
a penetration depth 5 ~ y/n/(uip) Stress effectively 
propagates instantaneously on length scales shorter than 
this. 

This picture generalizes to the case of homogenous 
viscoelastic media characterized by an isotropic, time- 
dependent shear modulus G(t) [23|], although the prop- 
agation of stress generally becomes super-diffusive [llj . 
We further assume that the medium is incompressible, 
which is a particularly good approximation for polymer 
solutions such as those considered here, at least at high 
frequencies [2(| [24|, [25| . The deformation of the medium 
is characterized by a local displacement field u(r , t) , and 
the viscoelastic analogue of the Navier-Stokes equation 
© is: 

P^u(r, t) = V • a(f, t) — VP + f(r, t) , (2) 

where 



II. THEORY AND BACKGROUND 

Newtonian liquids are described by the Navier-Stokes 
equation, which is non-linear. The non-linearity, how- 
ever, can usually be neglected either for small distances 
or for low velocities [l|, |2j. This is the so-called low 
Reynolds number regime, since the relative importance 
of non-linearities is characterized by the Reynolds num- 
ber Re = — where U, L, /?, and r\ are, respectively, 
the characteristic velocity and length scales, the density, 
and the viscosity. For steady flow, this regime can also 
be thought of as the non-inertial regime, in which stress 
propagates instantaneously and, for instance, the veloc- 
ity response at a distance r from a point force varies 
as l/r P, 0, [I[ . Such Stokes flow accurately describes 
the motion of micron-size objects in water on time scales 
longer than a few microseconds. 

Even at low Reynolds number, however, there are re- 
maining consequences of fluid inertia for non-stationary 
flows [3j. This unsteady Stokes approximation is de- 
scribed by the linearized Navier-Stokes equation: 



d _ 
P dt° 



rjW v — VP + /, 



(1) 



where v is the velocity field, P is the pressure that en- 
forces the incompressibility of the liquid and / is the force 
density applied to the fluid. By taking the curl of this 
equation we observe that the vorticity Cl = V x v satisfies 
the diffusion equation with diffusion constant v = n/ p. 

As noted above, the short-time response of a liquid 
to a point force generates a vortex. The propagation of 
stress away from the point disturbance is diffusive: after 
a time t, this vortex expands away from the point force 
to a size of order 5 ~ yfnt/ p. In the wake of this moving 
vortex is the usual Stokes flow that corresponds to a l/r 



a{r,t) = 2 / dt'G(t-t')j(r,t') (3) 



is the local stress tensor and 
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W+ I Vu 



(4) 



is the local deformation tensor. Incompressibility corre- 
sponds to the constraint V • u = 

Equations (12IMI) can be simplified by a decomposition 
of the force density and deformation into Fourier com- 
ponents. Taking spatio-temporal Fourier Transforms de- 
fined as 

u(k, u) = f d 3 r f dt e^-^u^t), (5) 

J J — OO 

and denning the complex modulus 

poo 

G*(oj) = G'(uj)+iG"(ui) = / dte luJt G{t), (6) 



we can eliminate the pressure by imposing incompress- 
ibility in Eqs. (03]). This leads to 



u(k,to) 



kk 



G*(cu)k 2 - ptu 2 



(7) 



where k — k/\k\. We invert this Fourier transform to 
obtain the displacement response function due to a point 
force applied at the origin. 



A. Response functions 



For a point force / at the origin, the linear response 
of the medium at r is given by a tensor otij, where 



3 



Ui(r,io) = aij(r,uj)fj(0,uj). Here, aij = a' i3 + ia" 3 is, in 
general, complex. Given our assumptions of rotational 
and translational symmetry, there are only two distinct 
components of the response. These are (1) a parallel re- 
sponse that is given by a displacement field u parallel 
to both / and r, and (2) a perpendicular response given 
by u parallel to / and perpendicular to r. The parallel 
response function aii, for instance, is obtained from the 
inverse Fourier transform of Eq. ([7]) [ll| . 

The response functions for general G*(u)) are given by 



a\\ (r, lu) = a'u + i ajj 1 



1 



47rG*(w)r 



X|| (ry/K) , (8) 



and 



a± (r, lo) = a ± + i a" = 



1 



87rG*(w)r 



XX (r^) , (9) 



where n = pui 2 /G*(uj) is complex and 



X\\ (x) 



[(1 - ix) e lx - 1] 



and 



X i_ (x) = — [l + (x 2 - 1 + ix) e 1 ' 



(10) 



(11) 



The magnitude of k defines the inverse (viscoelastic) pen- 
etration depth S. 



1. Simple liquids 

For a simple liquid, G*(lo) = —iuj-q and k = ipuj/rj. 
The velocity response of the liquid is then characterized 
by 



— iu>a\\ (r, to) = coa'u — i waj = 
and 

— iua±_ (r, u) — Lua'[^ — i Loa'j_ 
where 

xf (») 
X'± (x) 
X'i(x) 




[(1 + x) sinx — xcosa;] e x 
1 — [(1 + x) cos a; + x sin x] e" 



[(a; + 2a; 2 ) cosx — (1 + .x) 



sinx e 



(14) 
(15) 
(16) 



\[x + 2x 2 ) sinx + (1 + x) cos x] e x — l 



(17) 



Thus, for instance, the in-phase and out-of-phase velocity 
response in the parallel direction are given by ua" = 
-t— — y m and —loo 1 = — !— v{f. 
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FIG. 1: Comparison of theoretical and experimentally mea- 
sured response functions for a simple liquid. The predictions 
of the normalized velocity field from the dynamic Oseen tensor 
in Eqs. (|12I13|I are shown as black and gray lines. Normal- 
ized complex inter-particle response functions between two 
probe particles (silica beads, R = 0.58 (im) measured with 
the active microrheology method in water: (a) Anrrjuja\\ in 
the parallel direction and (b) 8nrr]LJct± in the perpendicular 
direction, plotted versus the ratio of the separation distance 
r (fixed for a given bead pair, r = 11.3 /im for parallel, and 
r = 10.3 /im for perpendicular) to the frequency-dependent 
viscous penetration depth 5 V . Both real (filled symbols) and 
imaginary parts (open symbols) are shown, for both paral- 
lel and perpendicular directions. These results are compared 
with the theory for r\ — 0.97 mPa s and p = 1000 kg m~ 3 . 
There is good agreement with no free parameters. 



We have written the response functions in Eqs. ([8j[9]) in 
a form in which the noninertial limits [x — > 0) are simple: 
Xii j_ - * 1- Thus, for a simple liquid in the limit x — ► 0, 
Eqs. (fT2l fil{|) reduce to the (time-independent) Oseen 
tensor [2|, |4|. For finite x, these response functions give 
the dynamic Oseen tensor 0, [26| , which are shown as the 
solid lines in Fig. [2 where for small r/S the parallel and 
perpendicular velocity response (i.e., —iu>a\\±) approach 
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and for a unit force at the origin. These then 
decay for r > 5. The region of negative response in the 
perpendicular case corresponds to the back-flow of the 
vortex. 

The response functions above represent the ensemble- 
average displacements due to forces acting in the 
medium. These response functions also govern the equi- 
librium thermal fluctuations and the correlated fluctua- 
tions from point to point within the medium. The rela- 
tionship between thermal fluctuations and response is de- 
scribed by the Fluctuation-Dissipation theorem. Specif- 
ically, for points separated by a distance r along the x 
direction, 



(18) 



where 



and 



C»(r,u)= dte^{u x (S),0)u x {r,t)) (19) 



which has both solvent and polymer contributions. For 
the Rouse model of flexible polymers z = 1/2 [27J , while 
for semifiexible polymers z = 3/4 [2(il2gl,|29j|. The latter 
case is shown as the solid lines in Fig. [5] We see that 
the oscillatory or anti-correlated response becomes more 
pronounced in viscoelastic materials. 

The magnitude of the complex modulus is given by 

\G\ = \J (gu z f + (wry) 2 + 2cu z+1 r 1 g sin (ttz/2), (22) 
while its phase is given by 

(lot) + guj z sin (ttz/2)) 



and 



COS'0 



I CI 



goj z cos (ttz/2) 



(23) 



\G\ 



(24) 



C ± (r,u)= dte iut (u y (0,0)uy(r,t)). (20) 



It is also useful to have the following expressions for the 
half- phase-angles 



2. Polymer solutions 

An experimentally pertinent illustration is given by the 
high frequency complex shear modulus of a polymer so- 
lution, 



ip if gu> z cos (ttz/2) 

sm 2 = Vn 1 — m — 

V> 1 1 ( gu z cos (ttz/2) 

cos 2 = Vn 1+ — m — 



(25) 
(26) 



G*(uj) = -iut) + g(-iuj) z = \G\e 



-iip 



(21) 



We define the real parameter (3 = ry pu 2 /\G\ and use the definitions, Eqns. (jlOllip to obtain the following compact 
expressions for the response functions (aj_,a||) which can be expanded using the compound angle formulae and the 
definitions above. The real and imaginary parts of the parallel response are given by 



4:Tr\G\raUr,uj) 



J 2 



and 



4tt 



|G|rof(r )W ) = -^ 



-j3 sin - 



-j3 sin ■ 



<-<^ { f3 cos — ) - ■ > ! sili 



sin j3 cos 



(3 cos 



^ , a ^ 
hp cos — 

2 A 2 



^ , a ^ 

hp COS — 

2 A 2 



while the corresponding expressions for the perpendicular response are given by 



%TT\G\ra' ± (r,u) = — \l 



-(3 sin ^ 



cos j3 cos — J + j3 sin I — + cos — J — j3 cos [ip + [3 cos — 



and 



%n\G\ra'[_(r,w) = -^e^ 8 ' 11 ^ |-sin f(3cos^\ +/3cos ^ +(3cos^j + (3 2 sin ^ + /3cos| 



(27) 



(28) 



(29) 



(30) 



The imaginary part of the response functions will be used to calculate the correlation functions, Eq. (|18p . used for 
analysis of the passive experiments whilst the real part of the response functions will be used for comparison with the 
active experiments. 
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We can simplify the expressions above in the limit that the polymer contribution to the viscoelasticity dominates 
the shear modulus. We then obtain the simple scaling form G*(oj) ~ g(— «w) z [30J giving \G\ = g\uj\ z ,ip — ttz/2. 
Further simplification of Eq. (|18p using the expressions in Eqs. (|3QI28p and definitions in Eqs. (|10|lip leads, e.g., to 



C,|(r,w) 



2T7u\G\r I/? 2 



{ 4 e " Sin( ^ )/? K 1 + sin (x) P) sin [ cos (x) 0] " cos (f) /3 cos t cos (-) /?] ] }» ( 31 ) 



where /3 = puS 1 j\G\ characterizes the overall decay of stress due to inertia. This decay corresponds to super- 
diffusive propagation of stress for viscoelastic media with G ~ u> z and z < 1, since the response is limited to a spatial 
range that grows with time as t( 2 ~ z )/ 2 . 



The resulting displacement field, exhibiting the vortex 
pattern, is shown in Fig. [3] for a point force at the ori- 
gin pointed along the re-axis. This flow pattern exhibits 
specific inversion symmetries: v x (v y ) is symmetric (anti- 
symmetric) for either x — > — x or y — > —y, as can be seen 
by the fact that the (linear) response must everywhere 
reverse if the direction of the force is reversed. 



B. Velocity autocorrelations and the long time tail 

The self-sustaining back-flow represented in Fig. [3] 
gives rise to long-lived correlations that, for instance, af- 
fect the crossover from ballistic to diffusive motion of a 
particle in a liquid. For a simple liquid, the fluid velocity 
(auto) correlations {v(0, t) ■ v(0, 0)) decay proportional to 
~ |i|~ 3 / 2 . This is known as the long time tail [H, HI, HI . 
For a viscoelastic fluid, stress propagation is faster than 
diffusive, resulting in a more rapid decay of velocity cor- 
relations. The decay is, however, still algebraic. 



1. Simple liquids 
For a simple liquid, Eq. fl} means that 



-ibjp - 



■qk 2 ) Vi = U 



kikj f. 



(32) 



for the Fourier transforms. This gives the response in 
velocity Uj to a force component fj that is (thermody- 
namically) conjugate to a displacement Uj. We denote 
this response function by Xviu^ where 



Sij hi kj ^ 



(33) 



r/k 2 — iup 

The fluctuation-dissipation theorem then tells us that 
1 d 



kTdt 



(«i(f, 4)^(0,0)), (34) 



where this is valid only for t > because of causality 
in the response. The correlation function is, however, 
defined for both positive and negative times. 

Due to translation invariance in time, the correlation 
function (vi(r,t + t')uj(0,t')) must be independent oft'. 



Thus, 



= ^(Vi(f> t + t') Uj (0,O) 



(35) 



= (vi(r,t + t') Uj (0,t')) + (vii^t + t^Vj&t')), 

which also means that 

kT Xvi u 3 {r,t)) = ( Vl (f,t)v,(0,0)). 

Again, this is valid only for t > because of causality in 
the response. Ultimately, however, we are interested in 
the autocorrelation function (v ■ v) = (viVi), for r — * 
above. In this case, the correlation function is manifestly 
symmetric in t. Thus, 



(vi(r -» 0,^(0,0)) = kT Xvt u t (r^ 0, \t\). 



(36) 



We first calculate XviUj (k, tj , since the limit 
XviUi (r— > 0,i) can be obtained from this simply by in- 
tegrating over all k. 



Xv t u J (k,t) = (Sij-kikj^j J — < 

\Sij k{kj 

—lixip 



I 



rjk 2 — iuip 



lu + ivk 2 



where v = r//p. But, the last integral can only depend 
on the combination tvk 2 , since we can replace u> by C,vk 2 , 
where £ is dimensionless. Specifically, 



XviUj (k,t) 



{^ij kikj 



—2nip 

ij kikj 

p 



-iC,vk 2 t 



C + i 



-vk't 



(37) 



for t > 0. Otherwise, the result is zero. This integral can 
be done by integration along a closed contour containing 
the real line in either the upper half-plane for t < 0, or 
the lower half plane for t > 0. 
Finally, to get the limit 



XviuAr-^ 0,t) 



(38) 
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FIG. 3: The displacement field displays a clear vortex-like 
structure. Here, a force in the x direction is applied at the 
origin (as shown by the filled circle and arrow). Distances 
are shown in units of the penetration depth 8 — \/\G\/ '(pco 2 ). 
This example has been calculated for the Rouse model with 
z = l/2. 



Again, this is only for t > 0. Thus, 

(v(r^0,t)-v(0,0)) = — UTrv\t\r 3/2 . (40) 

P 



FIG. 2: Comparison to Eqs. (|8llll[) of both real and imagi- 
nary parts of the normalized inter-particle response functions 
between two probe particles (silica beads, R = 0.58 pm), mea- 
sured with the active microrheology method for two separa- 
tion distances r in 1 mg/ml entangled F-actin solutions, (a) 
47rr|G|a||(u>) and (b) 8nr\G\a±(u}), plotted versus the ratio 
of the separation distance r (r — 11.3 pm in parallel and 
r — 10.3 /jm in perpendicular direction) to the frequency- 
dependent viscoelastic penetration depth 8 ve . Here, the pa- 
rameters g and z in Eq. (|2ip were varied to obtain simulta- 
neous fits of all data sets to Eqs. (|27H30[) . using a single set of 
parameters z and g, while accounting for the solvent (water) 
viscosity. Both the real (filled symbols) and imaginary (open 
symbols) parts of both parallel and perpendicular response 
functions are in a good agreement with Eqs. (|27I30|I with op- 
timal parameters g = 0.22 ± 0.05 Pa s z and z = 0.78 ± 0.02. 
The corresponding theoretical lines are shown for z — 0.75 
and g = 0.22. 



we simply integrate: 



d 3 k 



/k 2 t 



(2tt)3 
(Ain't) 



-3/2 



(39) 



2. Viscoelastic media 

For viscoelastic media, the calculation is similar, ex- 
cept that 



Xv iUj (k,uj) = -iuJXu tUj (k,uj) 



Xvu(k,uj) \ 5ij - hkj j . 



where Xuiuj — and where we have defined 



Xvu{k,u) = - 



p \u 2 - G*k 2 /p 



(41) 



(42) 



for simplicity. As above, all the singularities in this must 
lie in the lower half-plane in order for the response func- 
tion to be causal. 

Evaluation of the inverse Fourier transform of Eq. (f4"2"l) 
to obtain Xvu(k,t) can be done by use of the Mittag- 
Lcfflcr functions. The Mittag-Lefner functions [3ll [32j. 
E a (z), which are entire functions parameterized by a con- 
tinuous parameter a can be defined by the power series 



E, 



fc=0 



(1 + ak) 



a > 



(43) 
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Note Ei(x) = exp(a;). Straightforward manipulation [3l[ 
of the definition, Eq. (|4"3")l show that their causal Fourier 
transforms are given by 



dte luJt O(t)E a (-at a ) 



luJ 



lo 2 — a (— mj)" 



- (44) 



where 6(i) is the Heaviside step function. Performing 
the inverse Fourier transform, we obtain 

Xvu(k,t) = E 2 ^ z [-(gk 2 /p)t 2 ~ z ] ; <>0 (45) 

for G* — g(—iuj) z . The asymptotic expansions for the 
Mittag-Lefner functions [3l|, US are 



1 



E a (-z) 



T(l + a) 



+ 0(z 2 ), z<l 



l (46) 
+ 0(z- 2 ), z>l 



T(l - a 

Thus, the velocity correlation function is given by 



(«(r->0,t)-«(0,0)> = 
2kT 4tt r 1/a 



(47) 



P (^) 3 Jo 



k 2 dkE 2 _ z \-{gk 2 /p)\t\ 



2-z] 



From the asymptotic properties of E a (t) , it is clear that 
the integral does not converge at k — > oo necessitating a 
finite cut-off (which we choose as the size of the probe 
particle). We can express the correlation function as 

, ^, m , 3 in „\ 2kT ( p\ 3 ^ 2 2 — z 



P W (2^) 2 



and A a = (\t\g/pa 2 y/ a . 

We can get an approximation to the value of I a by 
splitting the integral into two sections 



1 pA a 



I a = J dxE<<(-x a ) + J dxE^(-x 

using the two asymptotic forms for E a (t). 
We then obtain finally the expression 

(v(f -> 0,t) • v(0,0)) ~ Ci|t|- 3(2 - z)/2 + C 2 |t|- (5 - 3x)/2 

(49) 

where 



Ci 



4fcT /p 



(27T)2p 



3/2 



(50) 



and 



5r(3 - z)Y{z - 1) - 3r(z - 1) - 15r(3 - z) 
15r(3 - z)T[z - 1) 

4fcT /p\ 1 



C 2 = 



(27r) 2 ap Vff / r ( z - !) 



(51) 



Here, the dominant first term in Eq. (|49[) corresponds to 
a faster asymptotic decay of correlations for z < 1 than 
for simple liquids. This is a direct consequence of the 
super-diffusive propagation of stress in this case. Also, 
it is interesting to note that the second term in Eq. ([49]) 
strictly vanishes in the z — > 1 limit of simple liquids. 



III. MATERIALS AND METHODS 

Simple liquids: For our experiments we used two 
Newtonian fluids with different viscosities r\ and mass 
densities p, namely water (77 = 0.969 mPas and p — 
1000 kgm -3 ) and a (1 : 1 v/v) water /glycerol mixture 
(77 = 6.9 mPas and p = 1150 kgm -3 ). 

Viscoelastic fluids: We performed experiments with 
two different viscoelastic fluids, worm-like micelle so- 
lutions and solutions of the cytoskeletal biopolymer 
F-actin. Worm-like micelles were prepared by self- 
assembly of cctylpyridinium chloride (CPyCl) in brine 
(0.5M NaCl) with sodium salicylate (NaSal) as countcri- 
ons, with a molar ratio of Sal/CPy = 0.5. Three differ- 
ent concentrations of worm-like micelles were used: c m = 
0.5%,1% and 2%. Worm-like micelles behave essentially 
like linear flexible polymers with an average diameter of 
about 3 nm, a persistence length of about 10 nm, and 
contour lengths of several /_tm [34| . F-actin was polymer- 
ized from monomeric actin (G-actin) isolated from rab- 
bit skeletal muscle according to a standard recipe [35| . 
G-actin was mixed with silica beads in polymerization 
buffer [2 mM HEPES, 2 mM MgCl 2 , 50 mM KC1, 1 mM 
Na 2 ATPa, and 1 mM EGTA, pH 7] and incubated for 
1 hour. Entangled actin solutions were used as a model 
system for semifiexible polymer solutions. Experiments 
were done at concentrations of c = 0.5 and 1 mg/ml. 



A. Experimental methods 

Details of the experimental set-up can be found in 
[HI, HH; HE H3 ■ Briefly, we have used a custom-built in- 
verted microscope that includes a pair of optical traps 
formed by two focused laser beams of different wave- 
lengths (A = 1064 nm, ND:YV04, Compass, Coherent, 
Santa Clara, CA, USA) and A = 830 nm (diode laser, 
CW, IQ1C140, Laser 2000). The optical traps fulfil two 
functions: (i) They confine the particles around two well- 
defined positions and at the same time detect particle 
displacements with high temporal and spatial resolution 
in the passive method, (ii) In the active method one trap 
is used to apply a sinusoidally varying force to one par- 
ticle, while the response of the other particle is detected 
in the second trap. In the passive method, a pair of 
silica beads (Van't Hoff Laboratory, Utrecht University, 
Utrecht, Netherlands) of various radii (R = 0.25 pm 
±5%, 0.580 pm ±5%, 1.28/jm ±5% and 2.5 pm±5%) 
were weakly trapped (trap stiffness between 2 pN/m and 
5 /_iN/m, where larger particles required the higher laser 
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intensities to avoid shot noise. Transmitted laser light 
was imaged onto two quadrant photo diodes, such that 
particle displacements and in the x and y directions were 
detected interferometrically [12 ]. A specialized silicon 
PIN photodiode (YAG444-4A, Perkin Elmer, Vaudreuil, 
Canada), operated at a reverse bias of 110 V, was used 
in order to extend the frequency range up to 100 kHz for 
the 1064 nm laser [38(. The 830 nm laser was detected 
by a standard silicon PIN photodiode operated at a re- 
verse bias of 15 V (Spot9-DMI, UDT, Hawthorne, CA). 
Amplified outputs were digitized at 195 kHz (A/D inter- 
face specs ) and further processed in Labview (National 
Instruments, Austin, TX, USA). Output voltages were 
converted to actual displacements using Lorentzian fits 
to power spectral densities (PSD) as described in [39|. In 
the case of water, calibration was done on the beads that 
were used in the experiments, while for the viscoelastic 
solutions and the more viscous liquid, calibrations were 
done in water with beads from the same batch. 

In the active method, the 1064 nm laser was used 
to oscillate one particle, while the 830 nm laser was 
used for detection of the second particle at a sep- 
aration distance r. The driving laser was deflected 
through an Acousto-Optical Deflector (AOD) (Te02, 
Model DTD 276HB6, IntraAction, Bellwood, Illinois), 
using a voltage-controlled oscillator (VCO) (DRF.40, AA 
OPTO-ELECTRONIC, Orsay, France). The force ap- 
plied to the driven particle was calibrated by measuring 
the PSD of the Brownian motion of a particle of the same 
size trapped in water with the same laser power [3!|. 
The output signal from the QPD detecting the second 
laser was fed into a lock-in amplifier (SR830, Stanford 
Research Systems, Sunnyvale, CA, USA) to obtain am- 
plitude and phase of particle response. All experiments 
were done in sample chambers made from a glass slide 
and a cover slip with about 140 /im inner height, with the 
particles at at least 25 //m distance from both surfaces. 
The lab temperature was stabilized at T = 21.5 °C. 

IV. DATA ANALYSIS 

In both the active and passive methods we calculate 
the linear complex response function a defined by u(to) = 
a(uj) x F(u), where F(u>) is the applied force. Linear re- 
sponse applies by definition in the passive method, and in 
the active method the particle displacements u{ui) were 
kept sufficiently small. Again, we consider separately 
real, a' (to), and imaginary parts, a"(u>), of the response. 
In all our experiments, as sketched in Fig. |4] the coordi- 
nate system was chosen in such a way that x is parallel 
to the line connecting the centers of the two particles (||) 
and y perpendicular (-L)to that. The inter-particle re- 
sponse functions along these two directions were used to 
determine the flow field. 

The displacement uP(u>) of particle 1 in the x di- 
rection is related to the force Fx acting on particle 
2 according to (a;) — a\\(w) x F^\uj). Similarly 




FIG. 4: (Color online) Schematic sketch of the experiment. 
A pair of silica beads (radius R) is trapped by a pair of laser 
traps at a separation distance r. In the passive method, the 
position fluctuations of each particle in the x and y directions 
are simultaneously detected with quadrant photodiodes and 
the displacement cross-correlations are measured parallel and 
perpendicular to the line connecting the centers of the two 
beads. In the active method, one of the beads is oscillated 
by rapidly moving one laser trap either in x or in y direction 
and the resulting motion of the other particle is measured in 
x and in y direction. Laser intensity was adjusted to result in 
the trap stiffness of typically between {2y,N/m and 5/j.N/m) 
for passive measurements. 



the perpendicular response function was derived from 
Uy(u)) — aj_(uj) x Fy 2 \uj). The single-particle re- 
sponse functions for each x and y directions are de- 
fined as Ux^(u) — a au to(w) x Fx(cj). For homoge- 
neous, isotropic media, the two functions aii com- 
pletely characterize the linear response at any point in 
the medium due to a force at another point. The dis- 
placement response functions a||.j_(u;) determine both 
position and velocity response —iu>a\\±(u>). 

In the passive approach, the medium fluctuates in 
equilibrium, and the only forces on the particles are 
thermal/Brownian forces. Therefore the fluctuation- 
dissipation theorem (FDT) of statistical mechanics [4(| 
relates the response of the medium to the displacement 
correlation functions. For two particles, these correla- 
tion functions are the cross-correlated displacement fluc- 
tuations: (u^tyu^iO)) and (uPtyu^tp)}. We used 
Fast Fourier Transforms (FFT) to calculate displacement 
cross-correlation functions in frequency space and ob- 
tained the imaginary parts of the complex inter-particle 
response functions a!,', j_(uj) via the FDT: 
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where k is the Boltzman constant and T is the controlled 
laboratory temperature. The real parts of inter-particle 
response functions aj, were obtained by a Kramers- 
Kronig integral: 



<4xM 



-p 



(a" 



I.-L 



( W ) 



c 2 - 

cos(iw) 



(54) 



(C) sin(*CR, 



where P denotes a principal- value integral [20] . The high 
frequency cut-off of the Kramers-Kroning integral lim- 
its the frequency range of the calculated a',, ± (u>) [HI]. 
We also used the active method to obtain both real and 
imaginary parts of the response functions with 100 kHz 
bandwidth, as in Refs. [U, [3?j]- Here the lock-in am- 
plifier provides directly in-phase (real part) and out-of- 
phase (imaginary part) response of the second particle. 
The measurements were done over a grid of driving fre- 
quencies. 



V. RESULTS 

A. Simple liquids 

In the low- frequency limit, where fluid inertia can be 
neglected, the inter-particle response functions are in- 
versely related to the shear modulus of the medium 
[20l I24I l4lj . For a simple viscous fluid, the response 
functions in this limit are given by: 
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where r is the separation distance between the two par- 
ticles and r\ is the viscosity. These relations also shows 
(via the FDT) a 1/r dependence of the Fourier transform 
of the displacement cross-correlation functions: 



and 
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In Fig. [5l displacement cross-correlation functions for 
particle pairs in water, normalized to compensate for the 
distance dependence (4irrS\\ and 8irrS± ), are plotted 
versus frequency. For comparison a single-particle dis- 
placement auto-correlation function, normalized for the 
bead size dependence (6irRS aut o), is plotted as a solid 
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FIG. 5: Normalized displacement cross-correlation functions 
(a) 47iTiS'|| and (b) 8nrS± of two probe particles (silica beads, 
R = 0.58 fim) in water versus frequency (/ = lj/2h), com- 
pared for different separation distances r (symbols). In both 
(a) and (b), the solid line is the auto-correlation function of 
a single particle, normalized by 6nR. The single-particle mo- 
tion agrees well with the expected frequency dependence of 
a single Brownian particle (slope of -2), as indicated by the 
dashed lines. 



line, for a particle radius of R = 0.58 fim. The auto- 
correlation function agrees well with the power-law slope 
of -2 up to nearly 100 kHz. This slope is expected 
from the high-frequency limit of the Lorentzian shape of 
the power spectral density (= Fourier Transform of the 
displacement autocorrelation function) of the displace- 
ments of a harmonically confined Brownian particle in 
a viscous fluid 39]. The effect of fluid inertia is evi- 
dent as a deviation from this power law in the displace- 
ment cross-correlation functions. A systematically r de- 
pendent decrease of the cross-correlations is apparent at 
high frequencies for separations r ranging from 2.2 /im 
to 11.7 /im. The faster decrease of cross-correlations is a 
manifestation of the finite velocity at which stress prop- 
agates into the medium. For larger r, the data show 
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that the decrease begins at a lower frequency, because 
it takes longer for stress to propagate further. A com- 
parison of Figs. and b shows that the decrease of 
the cross-correlation is, at the same separation distance, 
more pronounced in the perpendicular channel than in 
the parallel channel. This is due to the fact that in the 
vortex-like flow pattern of Fig. [3] there is a region of fluid 
motion in the opposite direction to the applied force. For 
r > 5 fim (open squares), the cross-correlations become 
negative in the observed frequency window (not shown 
in the log- log plot). At still higher frequencies, the dis- 
placement cross-correlation functions again become pos- 
itive (Fig.[5jD), which is visible for the larger separations, 
r = 8.3 fim and r = 11.7 fim, consistent with the ex- 
pected oscillation in the displacement cross-correlation 
functions in the frequency domain. This effect becomes 
more pronounced in viscoelastic media 

The spatial and temporal propagation of the inertial 
vortex in a general viscoelastic medium is characterized 
by Eqs. (|gl[TTj) 0,0, EH. For simple liquids, A = rj (vis- 
cosity) and z = 1. In Figs. [5^ and b, we compare the 
imaginary parts of the normalized inter-particle response 
functions <k]rqru)Ci'n {uj) (Fig. [5^) and &irrjruja" L {uj) (Fig. 
[BJd) for water and a (1:1 v/v) water/glycerol mixture. In 
order to collapse all data onto a single master curve, as 
suggested by Eqs. (|8l fTT|) . we have plotted these normal- 
ized response functions versus the probe particle separa- 
tion r scaled by the corresponding frequency-dependent 
penetration depth 5 V = \J rffpuj. As shown in Fig. [SI data 
taken at all of the different separations r ranging from 
2.2 fim to 11.7 fim fall onto a single curve for both, the 
parallel (Fig. [6^) and the perpendicular (Fig. [6Jd) inter- 
particle response functions. Data for water also collapse 
on water/glycerol data, after accounting for the different 
viscosities, which are known in both cases. Thus, no free 
fit parameters were used. The single curves of collapsed 
data are in quantitative agreement with the frequency- 
dependent dynamic Oseen tensor in Eqs. (|12I13|) shown 
by the solid lines in Figs. [Hk and b, where the region of 
negative response in the perpendicular direction corre- 
sponds again to the back- flow region of the vortex already 
described above. 

At separations r large compared to particle size, the 
inter-particle response functions become independent of 
probe particle size and shape [24j, leaving a dependence 
only on r and uj. In our analysis, we have assumed that 
the particles are point-like since the ratio r/R > 4 in all 
experiments. In order to directly check the validity of the 
approximation, we measured the inter-particle response 
functions with particles of different sizes. In Fig. [7^, and 
b the imaginary parts of the normalized response func- 
tions in water Airr]rujot'^{uj) and STirfruja'^ {uj) are plotted 
versus r/S v , obtained with probe particles of radius R = 
0.58 fim, 1.05 fim, 1.28 fim and 2.5 fim. We find no sys- 
tematic bead-size dependence, justifying the point-probe 
approximation. Again all normalized and scaled data col- 
lapse onto one single curve for each channel, which in turn 
agrees well with the dynamic Oseen tensor. The devia- 




FIG. 6: Normalized imaginary parts of inter-particle response 
functions between two probe particles (silica beads, R — 0.58 
fim) measured with the passive method, (a) inrijuia'^ in the 
parallel direction and (b) 8nrrjuja'/_ in the perpendicular di- 
rection, plotted versus the ratio of the separation distance 
r (fixed for a given bead pair) to the frequency-dependent 
viscous penetration depth 5 V , in water (open symbols, rj 
= 0.969 mPa s) and in water/glycerol (filled symbols, r\ = 
6.9 mPa s). Solid lines are Oseen's predictions for a simple 
liquid with no adjustable parameters. Data is only plotted 
for uj > 200 rad/s and for uj < 2krad/s one in every 5 data 
points is shown. 



tions observed in the last data points for R = 1.05 fim 
and 2.5 fim are most probably due to the influence of 
shot noise at high frequencies. Shot noise becomes a 
problem for larger beads because of their smaller fluctu- 
ation amplitudes, which will eventually produce signals 
that approach the fixed shot noise level. The deviations 
observed for these particle sizes are not consistent with 
a possible error due to finite particle size, which should 
be larger for smaller r/R. 

So far, we have considered the passive fluctuations that 
directly measure only the imaginary (out-of phase) part 
of the inter-particle response functions (ajj j_(oj)). With 
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FIG. 7: (Color online) Normalized imaginary parts of inter- 
particle displacement response functions between two probe 
particles of various radii measured with the passive method 
in water (silica beads, R given in legend), (a) 4nrrjuia"\ in the 
parallel direction and (b) 8nrgcja'^_ in the perpendicular di- 
rection, plotted versus the ratio of the separation distance r 
(fixed for a given bead pair) to the frequency-dependent vis- 
cous penetration depth 5 V . These response functions also rep- 
resent the in-phase velocity response normalized by the corre- 
sponding components of the Oseen tensor. Different particle 
sizes (R = 0.58, 1.05, 1.28 and 2.5 /im) were used at various 
separation distances r. 



the active method described in Sec. MI Al we can deter- 
mine both real and imaginary part of the response func- 
tions. In Figs. [1^ and b, we show the normalized inter- 
particle response functions a'^ and ajj ± (lo) mea- 
sured by active microrheology (bead radius R— 0.58 /im) 
in water for both parallel and perpendicular direction. 
In both cases we find good agreement with Eqs. (|12I13[) . 
The slightly different separation distances in parallel (r 
= 10.3 /im) and perpendicular (r = 11.3 fim) were due 
to different settings of the AOD signal in these measure- 
ments. 



B. Viscoelastic solutions 

In viscoelastic polymer solutions, the elastic compo- 
nent in the response of the medium modifies the propa- 
gation of the inertial vortex. We first discuss our results 
for worm-like micelle solutions, the viscoelastic proper- 
ties of which have been characterized by microrheology 
[42[ l43l| . In Fig.[8]we have plotted the displacement cross- 
correlation functions (AirrS^ and 8irrS±) in a 1% worm- 
like micelle solution versus frequency for two particles at 
various separation distances r between 2 and 8 \xm. For 
comparison, we have added the scaled auto-correlation 
function 6irRS au to for a single particle (R = 0.58 (im). 
In contrast to the situation in simple liquids, the particles 
were here confined by the surrounding polymer network 
and do not diffuse freely. Thus, the frequency dependence 
is weaker than for Brownian motion (i.e., the slope is less 
stee p th an -2 in the low-frequency, non-inertial regime) 
[42l |43|. As before, the displacement cross-correlation 
functions of the two probe particles are used to map the 
vortex-like flow pattern and its propagation in time. The 
r-dependent decrease of the cross-correlation functions 
occurs for both parallel (Fig.[S^,) and perpendicular (Fig. 
[8}d) directions, although the latter is more apparent. 

In order to collapse these data onto a single curve for 
each of the two channels (parallel and perpendicular), fol- 
lowing Eqs. (jSl fTTj) . we plot the normalized inter-particle 
response functions 47r|G|ra|j and 87r|G|m" versus scaled 
r/S ve , where S ve is the viscoelastic penetration depth. 
Unlike for the water and water/glycerol samples, we do 
not a priori know the frequency-dependent shear modu- 
lus G*{u>). Based on theoretical expectations for flexible 
polymers [27j , as well as on prior high-frequency rheology 
of worm-like micelle solutions (42l. |43T], we assume that 
the shear modulus has the functional form given in Eq. 
([2"Tj) . Thus, in order to achieve the collapse of all data 
onto the master curves represented by Eqs. (fSl fTTj) , we 
vary the two correlated parameters g and z. We expect 
z to be independent of the micelle concentration, while 
g should depend linearly on the polymer/micelle concen- 
tration. In Fig. [9^,, we show the resulting collapse of the 
normalized inter-particle response functions in parallel 
and perpendicular directions for different separations r. 
The predictions of Eqs. (|51 fTT|) are shown with black and 
gray lines ll|. The best overall collapse of the data for 
worm-like micelles solutions at all concentrations (0.5, 1 
and 2 weight percent) and separations r from 2 \xm to 
16 \xm was found for z = 0.68 ± 0.05 and g concentration 
dependence as seen in Fig. [SJs. 

In order to further test the inertial effects in viscoelas- 
tic media, we also performed experiments on another vis- 
coelastic fluid with a somewhat different frequency de- 
pendent shear modulus, namely entangled F-actin solu- 
tions. At high frequency, semiflexible F-actin filaments 
contribute to the viscoelasticity of the medium in a dif- 
ferent way from flexible polymers (l9l . [20L [28l . |29| . There- 
fore, the spatial structure and the propagation dynamics 
of the vortex should be different. Figure [TO] shows the 
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FIG. 8: Normalized displacement cross-correlation functions 
4nrS\ | (a) and 8irrS± (b) of two probe particles (silica beads, 
R — 0.58 /im) in worm-like micelle solutions (c m = 1 wt%) 
versus frequency (/ = lu/2tt) compared for different separa- 
tion distances r. The solid lines represent the auto-correlation 
function of a single particle normalized by 6irR. The dashed 
lines indicate slopes of -2, corresponding to diffusive motion. 



collapse of the inter-particle response functions 47r|G|raj| 
and 87r|G|ra" plotted versus scaled distance r/S ve onto 
two master curves for the parallel and the perpendicular 
direction. The actin concentration was lmg/ml and the 
probe radius 0.58 [im and we used separation distances 
r ranging from 4.2 fim to 16.2 fim. In an F-actin so- 
lution of this concentration, the magnitude of the shear 
modulus is large, therefore the vortex propagates faster, 
making it harder to observe. In particular, it was diffi- 
cult to determine the parameters g and z in this case. 
We found the best collapse with z — 0.78 ±0.1 for data 
taken in passive method. We then fixed z to 0.75 known 
from the power law dependence behavior reported previ- 
ously pi [H HI HI, and found g = 0.18 ± 0.13Pa s z . 
To reduce the large error bars in the passive method, it 
would be necessary to repeat our measurements at higher 
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FIG. 9: (Color online) (a) Collapse of the imaginary parts of 
the normalized inter-particle response functions (47rr|G|ajj (u) 
and 8nr\G\a'±(u>)) between two probe particles (silica beads, 
R — 0.58 /im), measured with the passive method for differ- 
ent separation distances r in worm-like micelle solutions of 
1 wt%, plotted versus the ratio of the separation distance r 
(fixed for a given bead pair) to the frequency-dependent vis- 
coelastic penetration depth S ve . Here, the viscoelastic pene- 
tration depths were determined by varying the parameters g 
and z in Eqs. (|8I11I21[ ) to obtain collapse. Optimal param- 
eters were z = 0.68 ± 0.05 and g = 0.0275 ± 0.008 Pa s z , 
where the solvent (water) viscosity has been taken into ac- 
count, (b) Dependence of (optimal) parameter g on micelle 
concentration for a fixed z. 



frequencies and/or larger separation distances. Neverthe- 
less, our results are consistent with prior measurements 
and predictions of both parameters. 

Independently, we have measured both real, aj, 
and imaginary, a'/, j_(u>), parts of the response functions 
directly by actively manipulating one particle and mea- 
suring the response of the other. In Figs. [2^i and b, the 
parallel (Fig. [2^l) and perpendicular (Fig. [2Jd) complex 
inter-particle response functions for c = 1 mg/ml F-actin 
solutions are shown, probed with beads of 1.28 fim ra- 
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FIG. 10: (Color online) Collapse of the imaginary parts of 
the normalized inter-particle response functions between two 
probe particles (silica beads, R = 0.58 /xm), measured with 
the passive method for different separation distances r in F- 
actin solutions of concentration 1 mg/ml, plotted versus the 
ratio of the separation distance r (fixed for a given bead pair) 
to the frequency-dependent viscoelastic penetration depth 
5 ve . Here, the parameters g and z in Eq. (|21|l were varied 
to obtain simultaneous collapse of all data sets onto Eqs. ((§]■ 
lll[) . using a single set of parameters z and g, while accounting 
for the solvent (water) viscosity. We find z — 0.78 ± 0.1, and 
g — 0.18 ± 0.13 as optimal parameters. Data are presented 
for parallel (closed symbols) and perpendicular (open sym- 
bols) directions. The solid black line represents Eq. ()28J | and 
the gray line represents Eq. (|30[1 . both with z = 0.75 and g = 
0.2 Pas 2 . 



dius. Here the inter-particle response functions were fit- 
ted with Eqs. (|8l [TTj) to find parameters g = 0.22±0.05 
Pa s z and z= 0.78±0.01 simultaneously. Data for both 
parallel (r = 12.1 /zm) and perpendicular (r — 13.5 /im) 
channels are compared with z — 0.75 and g — 0.22 Pa s z 
in Fig. [2^ and b. 

The values of z and g found from both methods are 
consistent and also agree with results from prior experi- 
mental microrhelogy and macrorheolog y exp eriments for 
both entangled actin solutions fl8l. [HA |20| . l28l [29j and 
worm-like-micelle solutions (42l. 1-431] - To obtain these 
values it was essential to model the inertial effects in- 
cluding both polymer and solvent contributions to the 
shear modulus. We observed that, although the high- 
frequency rheology of the polymer solution is dominated 
by the polymer, the background solvent contributes non- 
negligibly to the inertial vortex propagation. To test this, 
we excluded the solvent shear modulus (—iLurj) in Eqs. |[B} 
I11I21[) . and analyzed our data assuming a high-frequency 
shear modulus of the form G — guj z . For both worm- 
like micelle solutions and entangled F-actin solutions, we 
found much larger values of z ~ 0.9 and a nonlinear con- 
centration dependence of <?, contrary to expectations. 



VI. DISCUSSION 

In our experiments we have directly resolved the iner- 
tial response/flow of fluids on micrometer and microsec- 
ond time scales using optical trapping and interferometric 
particle-tracking. Our results demonstrate that vortic- 
ity and stress propagate diffusively in simple liquids and 
super-diffusively in viscoelastic media. One consequence 
of inertial vortex formation is the long-time tail effect ob- 
served in light scattering experiments f(| . To connect to 
these results, we calculated the velocity auto-correlation 
of a single particle from the displacement fluctuations. 
Unfortunately, the effect we are looking for is subtle and 
is difficult to detect in the presence of other factors. At 
the highest frequencies the vortex is still influenced by 
the finite probe size, and at intermediate frequencies the 
particle motion is already affected by the laser trap po- 
tential. The results were thus inconclusive. Similar prob- 
lems have been reported in Ref. @. The effects of in- 
ertia we have described here set a fundamental limit to 
the applicability of two-particle microrheology techniques 
which are based on the measurement of cross-correlated 
position fluctuations of particles [3, [H, [H, H3, EH ■ Iner- 
tia limits the range of stress propagation at high frequen- 
cies, stronger in soft media such as those studied here 
than in media with higher viscoelastic moduli. Inertia 
affects measurements at frequencies as low as 1 kHz for 
separations of order 10 micrometers, showing an appar- 
ent increase of the measured shear moduli below their 
actual values [43j]. Since the stress propagation is dif- 
fusive, or nearly so, even measurements at video rates 
can be affected for probe particle separations of order 
50 micrometers. As we have shown here, these inertial 
effects are more pronounced in the perpendicular inter- 
particle response functions than in the parallel ones. This 
suggests that one should obtain shear moduli from the 
parallel inter-particle response functions if one doesn't 
want to correct for inertia. In a more precise analysis 
of two-particle microrheology experiments, the fluid re- 
sponse function can not simply be modeled by a general- 
ized Stokes-Einstein relationship and has to be corrected 
for inertial effects according to the probed frequency as 
well as the particles separation. Such corrections, how- 
ever, will necessarily be limited, given the exponential 
attenuation of stress due to inertia. 
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